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Abstract. 

We have recently introduced the idea of making use of three-photon positron 
annihilations in positron emission tomography. In this paper the basic characteristics 
of the three-gamma imaging in PET are studied by means of Monte Carlo simulations 
and analytical computations. Two typical configurations of human and small animal 
scanners are considered. Three-photon imaging requires high energy resolution 
detectors. Parameters currently attainable by CdZnTe semiconductor detectors, the 
technology of choice for the future development of radiation imaging, are assumed. 
Spatial resolution is calculated as a function of detector energy resolution and size, 
position in the field of view, scanner size, and the energies of the three gamma 
annihilation photons. Possible ways to improve the spatial resolution obtained for 
nominal parameters: 1.5 cm and 3.2 mm FWHM for human and small animal scanners, 
respectively, are indicated. Counting rates of true and random three-photon events for 
typical human and small animal scanning configurations are assessed. A simple formula 
for minimum size of lesions detectable in the three-gamma based images is derived. 
Depending on the contrast and total number of registered counts, lesions of a few mm 
size for human and sub mm for small animal scanners can be detected. 
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1. Introduction 

Positron Emission Tomography (PET) is the method of choice in functional medical 
imaging, both in clinical practice and research involving small animals. It uses short- 
lived positron emitting radionuclides to mark biologically active substances which can be 
traced while being metabolised in the body by detecting co-linear annihilation photon 
pairs and reconstructing the image. 

Medical imaging is not the only application of positrons. Apart from basic nuclear 
and elementary particle research, the physics of positron annihilation investigates the 
interactions of positrons with matter (Charlton and Humberston 2001). They are a 
sensitive probe allowing to acquire information about subtle structural and chemical 
properties of materials by precise measurements of annihilation radiation. Three 
main experimental methods are employed: positron lifetime spectroscopy, angular 
correlation of annihilation radiation and Doppler broadening of annihilation radiation. 
Combinations of the above are also in use. None of these concepts is currently used 
in PET, apart from accounting for the non-colinearity of the annihilation photons and 
finite positron range contributing to deterioration of spatial resolution. Annihilation 
pairs are used merely to determine the activity distribution changing in time. 

Wc introduced recently the idea of three-photon PET imaging (Kacpcrski et al 
2004, Kacpcrski and Spyrou 2004a). It was shown that the relatively rare (about 0.5 % 
in water) 3-photon positron annihilations (Ore and Powell 1949, De Benedetti and Siegel 
1954) can also be used for imaging. By detecting the positions and energies of the three 
photons, one can easily locate the point where the annihilation occurred. Thus the 
amount of information obtained from a single event is higher than for a 27 pair, where 
localisation is along a line. The rate of 87 decay is not only proportional to concentration 
of activity but is also sensitive to the local physical and chemical environments, notably 
the presence of oxygen (Cooper et al 1967, Klobuchar and Karol 1980, Hopkins and 
Zerda 1990, Kakimoto et al 1990). This is due to formation of positronium: a positron 
- electron bound state, which behaves as an active chemical particle. 75 % of the 
positronium is formed as an ortho-positronium, triplet state which in vacuum annihilates 
only into 87 with a relatively long lifetime of 142 ns. In matter interaction processes, in 
general called quenching, that lead to 27 annihilations are usually much faster, reducing 
the fraction of three-photon annihilations. 

Measuring and imaging oxygenation of tissues is important in various clinical 
contexts (MachuUa 1999), e.g. in oncology (tumour hypoxia, radiosensitivity) , cardio- 
and cerebrovascular disease and infections by anaerobic microorganisms. Intensive 
research is ongoing to develop radiotracers for imaging hypoxia. When the 87 
annihilations, which are simply ignored in the current PET scanners, are recorded, 
the positron itself, or more precisely the positronium, could serve as an oxygen-sensitive 
tracer. Oxygen is known to be a strong positronium quencher, therefore hypoxic regions 
should be characterised by higher 87 rates than those well oxygenated. The oxygenation 
image would be obtained alongside e.g. a routine FDG PET, saving cost, patient 
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radiation dose and inconvenience incurred by extra specific scan. 

In this paper we explore further the idea of 87 imaging by assessing the dependence 
of spatial resolution on various system parameters hke the size of the detectors, their 
energy resolution, size of the scanner, position within the field of view as well as the 
particular combination of three photon energies. We focus on two scanner configurations 
commonly used in practice: a small animal and a full body human scanner. The expected 
87 counting rates are estimated and examples of images obtained from 87 annihilations 
are produced. A formula for the minimum size of detectable lesions with respect to 
contrast and number of counts recorded is derived. 



2. Theory 

Let us recall the basic principles of 87 PET imaging (Kacperski et al 2004). Consider 
a 87-decay event that occurs at a point r = (x, y, z) (see figure H}. 



Figure 1. Imaging by tliree-photon anniliilations. 

If the three annihilation photons have energies Ei, E2, E^ and are detected 
in coincidence by detectors at points ri,r2,r3, respectively, then, from the law of 
momentum conservation, we get 

El X — Xi E2 X — X2 E^ X — X3 

Px = — 1 r H 1 r H 1 r = 

c|r — ri| c|r — c|r — ral 

El y-yi E2 y-y2 E3 y - y^ 

Py = — 1 r + — 1 r + — 1 r = 

c|r — ri| c|r — r2| c|r — rsl 

El Z - Zi E2 Z - Z2 E3 Z - Z3 
Pz = 1 7 H 1 7 H 1 7 = 0. 

c|r — ri| c|r — c \r — 

(1) 
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In addition, the law of energy conservation has also to be satisfied: 

Ei + E2 + E3 = nieC^ ^ 1022keV. (2) 

The energies have a continuous spectrum from to 511 keV with the probability 
distribution approximately uniform in the E1E2 space (Ore and Powell 1949, Chang and 
Tang 1985). With known detector positions ri, r2, and r^, the measurement of photon 
energies Ei, E2, and E^, enables the solution of the nonlinear set of equations (PJ) to 
determine the point r at which annihilation took place. Since the energy resolution 
and size of the detectors are finite, the location of annihilation is broadened into a 
region surrounding the point r. Note, however, that in contrast to the 27 decay, we 
obtain (neglecting the finite measurement accuracy) full information on the position of 
annihilation from a single event, rather than just a line of response. Unlike the rather 
intricate image reconstruction methods required for 27 events, one only needs to solve 
the nonlinear equations to retrieve r from a 87 event. The image is thus formed as a 
set of dots, each corresponding to a single 87 positron decay. It can be processed further 
by performing appropriate attenuation correction and e.g. window averaging to obtain 
the usual grayscale or colour-coded pixels display. Let us stress that the recording of 
three-photon events is done simultaneously with the prevailing 511 keV photon pairs. 
87 imaging is by no means an alternative to conventional PET, but rather an add-on, 
making use of the annihilation radiation which is currently wasted, it provides extra 
information. Dividing the 87 image by the 27 one we obtain a map of 87/27 decay 
probability ratio. This is actually a new imaging modality. 

3. CdZnTe detectors for PET 

It is clear from the equation that uncertainty in the energies Ei, E2, and E^ 
will result in a spread of the reconstructed annihilation site r. As it is indicated in 
(Kacperski et al 2004) and section 14.1.11 of this paper, the energy resolution is actually 
the crucial factor determining the spatial resolution of 87 imaging. With currently 
dominating PET scanners which are based on scintillator detectors (energy resolution 
typically worse than 15%) it is not possible to obtain an accteptable 87 image. However, 
new detector materials are gradually making their way into PET imaging technology, 
notably the room temperature operating semiconductors CdZnTe (Moses et al 1994, 
Scheiber and Gaikos 2001, Verger et al 2004). The main reason for considering the 
new material is the improvement in intrinsic spatial resolution due to precise depth- 
of-interaction information, which results in better image resolution and is of primary 
importance, specially in small animal imaging (Stickel and Cherry 2005). Compact 
scanner design and much better energy resolution allowing efficient rejection of scattered 
events are further advantages. The latter factor also allows implementation of the 
3-photon modality. CdZnTe cameras are already becoming increasingly popular in 
SPECT. The price of the material, processing electronics and implementation of a new 
technology are still inhibiting factors, nevertheless, projects to built prototype small 
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animal CdZnTe PET scanners are already under way. Many of the material related 
problems, like slow charge collection, can be greatly reduced by the use of stacked thin 
position sensitive detector arrays (Nemirovsky et al 2001, Moss et al 2001, Redus et 
al 2004) with appropriate pulse processing electronics. Currently energy resolutions of 
1% at 662 keV can be achieved. Best timing resolutions range between 5 and 8 ns, and 
recently even 2.8 ns has been reported (Drezet et al 2004). 

4. Monte Carlo simulations: small animal and human scanner models 

Performance of PET imaging based on 3-photon annihilations has been investigated for 
two kinds of model scanners: small animal and whole body human. We assumed the 
usual cylindrical scanner geometry with square faced detector elements (no dead layers). 
The parameters of the scanners are given in Table H 

Table 1. Basic parameters of model scamiers. 

Parameter Small Animal Human 



Diameter 12 cm 80 cm 

FOV 15 cm 24 cm 

Diameter/FOV ratio 0.8 3.33 

Detector size 2 mm 4 mm 

Number of detectors in ring 188 628 

Number of rings 75 60 

Detector energy resolution at 662 keV 1% 

Minimum detected 37 energy {Emm) 150 keV 

Maximum detected 37 energy {Emax) 480 keV 



In order to investigate their influence on the spatial resolution chosen parameters 
have been modified during the simulations, with all others remaining fixed as in tabled 

Spatial resolution 

For each particular set of parameters at least 10^ detected three-photon events emitted 
from a point source placed in the centre of the scanner were simulated. This was 
assumed to obtain the point spreads with high statistical accuracy, and does not reflect 
the actually achievable counting rates of 87 photons; this issue is addressed separately 
in section 14.21 Finite size and energy resolution of detectors were the only sources of 
blurring in the reconstructed image. We neglected the range of positrons in matter and 
residual momentum of the electron-positron pair. These factors are briefly discussed 
in section El Any photon hitting the surface of a detector element in the scanner was 
assumed to be detected at the centre of that element. The maximum photon position 
error is therefore \/2/2 times the detector size. We assumed that the energy resolution 
of the detectors depends on photon energy according to 5E = A\/E, where the constant 
A was chosen so that the relative energy resolution at 662 keV is 1% (except section 
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I4.1.1|) . After the reconstruction procedure, as described in section |21 standard deviation 
and FWHM of the reconstructed point were calculated. 

4.I.I. Energy resolution of detectors Let us begin with the dependence of a point 
source blur on the energy resolution of the detectors which is shown in figure El Clearly, 



20 




detector energy resolution at 662 keV FWHM [%] 

Figure 2. Spatial resolution (spread of a point source) at the centre of the scanner 
as a function of the relative energy resolution of the detectors at 662 keV. 

good energy resolution of detectors is crucial for 3-photon imaging. With the scintillators 
currently used for which the energy resolution is usually worse than 15% at 662 keV 
one cannot obtain an acceptable image from the 87 events. However, with CdZnTe 
detectors of energy resolution 1%, or better, at 662 keV, the spatial resolutions that 
can be obtained become acceptable, although worse than typical values for conventional 
PET. At very high energy resolutions the influence of detector size can be noticed 
(saturation at a non-zero value for 6E = 0), whereas above 1%, energy resolution is the 
dominating source of blurring. 

4- 1-2. Position in the FOV In the following simulations the variation of the spatial 
resolution within the field of view has been investigated as the point source was moved 
along the radial and axial directions. The results for the radial direction are shown in 
figure El 

The point spread is a combination of errors in photon energy and position detection. 
In order to see the influence of both factors we ran separate series of simulations assuming 
perfect energy detection (error in photon position only), point-like detectors (error in 
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Figure 3. Spatial resolution (spread of a point source) as a function of the transaxial 
distance from the centre for the small animal and human scanners. 



energy only), and the realistic case when both errors are present. As one can see 
in figure El the energy error contribution increases, while the photon position error 
decreases with the transaxial distance from the central axis. In effect, the resolution is 
quite uniform across the FOV, deteriorating sharply only close to the scanner walls. For 
the human scanner the optimal resolution is actually achieved at about 25 cm off-centre. 
The variability along the z axis (figurelH) is even smaller, remaining below 2.5% and 
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4% for human and small animal scanner, respectively. The results of other simulations 
obtained for the point at the centres of scanners are therefore representative for the 
entire useful FOV. 




Figure 4. Spatial resolution (spread of a point source) as a function of the distance 
from the centre along the z axis of the scanners. 

It should also be noted that the shape of the point spread function is spatially 
nonuniform. For the point at the centre the spread in the transaxial direction equals 
1.06 cm and 0.2 cm FWHM, respectively, for human and small animal scanner, whilst 
in the z direction it is 0.26 cm and 0.136 cm, respectively. As we move off from the 
central axis the PSF becomes elongated in the radial direction. 

4.1.3. Size of detectors Keeping the geometric proportions of the scanners fixed as in 
tableHlthe size of (square) detectors was changed, and their number adjusted accordingly 
to keep the sizes of scanners constant. Figure El shows the results of the simulations 
with detector energy resolution 1% at 662 keV (solid lines) and perfect energy detection 
(dashed lines). 

The overall dependence, similar to figure El is close to linear. Comparison of both 
curves for each scanner reveals the contributions of the two sources of error and indicates 
the direction of potential improvements. For our model human scanner with detector size 
4 mm the reconstructed point spread results mainly from the energy detection error (cf. 
also figure ??), therefore decreasing the size of detectors, without improving their energy 
resolution, would hardly improve the spatial resolution. For the small animal scanner 
with 2 mm detector elements the contribution from photon position detection error is 
higher, so reducing effective detector size would bring some performance improvement, 
although energy resolution is still more critical. The influence of detector size becomes 
more important for points far from the centres of scanners (see figure ??). 
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Figure 5. Spatial resolution (spread of a point source) as a function of detector size. 
The dashed lines show the same dependance assuming perfect energy resolution (no 
energy blur). 

4.1.4- Size of scanner The influence of tlie scanner size on tlie spatial resolution was 
studied by changing the diameter D of the scanners and their axial length (number of 
rings) to preserve the geometric proportions as in tabled The size of detectors was kept 
constant and their number adjusted to cover the whole cylinder. Again, three separate 
simulations were run for each scanner: assuming perfect detection of photon energy, or 
position, and modelling errors in both. The results can be seen in figure IHl 

One can clearly see that the near-linear dependence is due to the energy resolution 
error which depends linearly on scanner size. The detector size component remains 
constant. For both human and small animal scanners reducing scanner size as much as 
possible is desirable from the point of view of 87 imaging. This dependence is the main 
reason why the spatial resolution of the small animal scanner is significantly better than 
that of human the scanner. 

4.1.5. Energies of photons and cut-off energies As it was mentioned in section |2l the 
three photons can have any combination of energies satisfying (j2)) and (^. The error of 
the annihilation site r recovered from is a rather complex function of Ei, E2 and E3. 
We have calculated the FWHM spreads of a point source at the centre of the human 
scanner emitting photons with a particular energy combination satisfying (j21) and ((Hi 
(figure IZj). The spreads were averaged over the remaining free parameters (directions of 
emission) . 

The resolution is optimal for the symmetric three-photon decay with all energies 
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Figure 6. Spatial resolution (spread of a point source) as a function of scanner size 
D for human (upper plot) and small animal (lower plot) scanner proportions. Solid 
lines: finite photon energy and position resolution, dashed: point size detectors, dotted: 
perfect energy resolution (no energy blur). 



equal 340.7 keV and deteriorates as we move towards extreme values. In practice we 
never detect energies from the entire 87 spectrum 0-511 keV. At the upper end the 
511 keV photopeak resulting from the dominating 27 decays has to be cut-off. On the 
other hand, there is a detector and noise related low energy detection limit. Figure [7| 
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Figure 7. Spatial resolution (spread of a point source) as a function of photon 
energies for the human scanner. 

indicates yet another reason to avoid registering 87 decays with extreme energies: they 
introduce large errors in position reconstruction. In figure |S] the point spread of a point 
source as a function of upper and lower 87 energy thresholds was plotted. The smaller 




Figure 8. Spatial resolution (spread of a point source) as a function of upper and 
lower energy threshold for 87 photons in human scanner. The dot indicates the energy 
window used in all simulations of this paper (cf. table ^ 

window around the optimal 340.7 keV we set, the better spatial resolution we achieve. 
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There is, however, the common sensitivity-resolution trade-off: as the window narrows, 
the number of counts detected drops, so a compromise choice has to be made. 

Figure also suggests a relatively easy way of improving spatial resolution without 
losing counts. In all our simulations we treat all the reconstructed 87 annihilation points 
equally. We could, however, using the relation shown in figure assign a weight to each 
point according to the particular energy combination, so that counts with energies close 
to optimal {Ei = E2 = E^ = 340.7 keV) would contribute more than those with energies 
close to extreme. Such approach would probably improve spatial resolution by about 
30 - 40 %. 

4-2. Counting rates and random scattered 37 events 

Poissonian noise is one of the key limiting factors in nuclear medicine imaging. With the 
probability of three-photon annihilations being two to three orders of magnitude lower 
than that of two-photon, the question of sufficient number of counts obviously needs to 
be addressed. Will it be enough to obtain meaningful images? 

We will assess the rates of 87 counts for our two model scanners with phantoms in 
the form of spheres filled with water and uniform positron activity distribution. The 
diameters of the spheres are 20 cm and 4 cm for the human and small animal scanners, 
respectively; they are placed at the centres of the scanners. The counting rate C^^ 
of true 37 decays is in general determined by probabilities of emission, detection and 
attenuation. It can be expressed by the formula 

Cs^ = A r3^/2^ Ecut V deff gsj (e"^" + giAre"^" + q2{AT)^e-^^), (3) 

where A is the total activity in the phantom. The meaning of the factors and their 
values for the two model systems is summarized in table El The expression in brackets 
approximates the probability that a true three-photon decay will not be obscured by 
(one or more) coincident 27 events. The three terms correspond to zero, one and 
two 27 pairs respectively. Higher order coincidences are negligible at the activities 
relevant for PET. The factors qi and q2 are the probabilities of unique identification 
of the true 87 event despite the coincident 27 photons. This is possible because the 
coincident photons may just miss the scanner, or, if they are detected, their energies 
are quite different from those of 87 annihilation. Numerical values of the parameters 
were obtained from simulations and analytical calculations. We assumed the scanners 
to be made of 3 cm thick CdZnTe detectors, which yields average full energy photopeak 
efficiency for the 87 photons of about 85%. The 87/27 ratio was derived from positron 
lifetime experimental data for water and organic liquids. In fact, it depends sensitively 
on the chemical composition of the solution. Usually dissolved ions as well as gases (in 
particular oxygen) lead to a decrease of 87 rate. On the other hand, the ratio increases 
in liquids of smaller surface tension. The r3^/27 ratio for pure, degassed water is about 
0.5%; in alcohols it reaches about 0.7%. The minimum value occurring in substances 
(e.g. metals) where no positronium is formed is 0.27%. Our choice of r3^/27 = 

0.4% 

seems therefore reasonable and not too optimistic. 
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Table 2. Parameters influencing 87 counting rate (cf. (jSJ). 



Symbol 


Parameter 


Small Animal Human 


V 

931 


attenuation factor 
geometry factor 


0.65 
0.51 


0.11 
0.052 


91 


probability that one 27 pair 

will not interfere with 37 detection 


0.94 


0.98 




probability that two 27 pairs 

will not interfere with 87 detection 


0.5 


0.8 


''37/27 


87/27 emission probabihty ratio 




0.004 


Ecut 


fraction of 87 photons in the 
[E„un, Emax) window 




0.655 


deff 


detection efficiency 




0.85^ = 0.614 



The 87 counting rates for both scanners are plotted in figures IHl and ^| They 
can amount up to about 1000 cps for human scanner and 80,000 cps for small animal 
scanner. The difference is mainly due to the larger solid angle covered and much smaller 
attenuation for the small object. 

In conventional 27 PET, apart from the true coincidence counts forming the image, 
there are scattered and randomly coincident photons contributing to image noise. In 
the case of 87 imaging those two kinds of events have to be considered simultaneously. 
Because of the resolution requirement we need to use high energy resolution detectors 
(see section HIlT}. The unique properties of 87 decay, in particular the energy condition 
(j21) enables then to efficiently distinguish the true events from accidental coincidences. 
Nevertheless, it can happen that all the conditions are fulfilled by chance leading to 
spurious counts. Let us estimate the rates of such counts for our model scanners. 
Because of the small 87 decay rate the probability of two or more 87 events occurring 
and being detected within the resolving time is very small. The main source of spurious 
counts are the randomly coincident 27 annihilation photons that are scattered or 
partially detected, so that they fall within the energy range of the 87 spectrum. Let 
us denote by P{E\511) the probability of a 511 keV photon emitted from the phantom 
to deposit the energy E in a. detector. It is actually a convolution of the spectrum 
of photons leaving the phantom (primary and scattered) and the detector response 
function. The probability of three coincident photons randomly fulfilling condition (j2I) 
within the resolution determined window AE can be approximately written as 

Ii= [ P(Ei\511)P(E2\511)P(E3\511)dEidE2dE3 (4) 

JwiAE) 

where the integration is over the region W{AE) between the planes Ei + E2 + E^ = 
1022 keV ± AE in the E1E2E3 space, further limited by Emin and Emax, the 
upper and lower hmits of detected 87 spectrum. Note that Ji does not include any 
geometry dependence. In our estimation P(£'|511) in the interval {Emm, Emax) can 
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be approximated by a linear function with parameters obtained from Monte Carlo 
simulations. The rate of false 87 counts can be in general expressed as 



The fraction is the Poissonian probability of n + 1 decays occurring within the resolving 
time r, and /(n) is the probabihty that at least one triplet of photons will be detected 
and falsely accepted as true 87 event. Practically, it is sufficient to consider just the 
first two terms in and then the false 87 counting rate becomes 



The two terms in the curly brackets correspond to two and three two-photon 
annihilations occurring within the resolving time, gi is the probability to hit the scanner 
by a single photon (geometric factor). We assumed that the directions of photons are 
uncorrelated, which is not quite true, but works well as an approximation. Even when 
three detected accidental photons fulfil condition (0) ±Ai?, when inserted into (Q) most 
of them would either give no real solution, or generate a point outside the object. Pobj 
is the fraction of events yielding a point within the object, contributing to 87 image 
noise. For our scanners and phantoms Pobj = 0.02, the single photon geometric factors 
5'i=0.287 for the human and ^'1=0.781 for the small animal scanner. 
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Figure 9. Estimated true and random 87 counting rates for liuman scanner. 
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Figure 10. Estimated true and random 87 counting rates for small animal scanner. 

The rates of random 87 events f figures lUl and fTUj) remain very low compared to true 
87 even at high activities, so they do not contribute significantly to image noise. This 
signal-to-noise ratio is, however, sensitive to the detection efficiency of the detectors and 
may worsen significantly for photopeak efficiencies lower than those assumed. 

4.3. Examples of images and detectahility of lesions 

In order to generate examples of images we used the human scanner with the phantom 
described in the previous section and added lesions in the form of spheres of variable 
size, lined up across the centre of the phantom as shown in figure fTTl fa). We define the 
lesion contrast as c = (ri/ — nh)/nh, where ni is the density of registered 87 counts in 
the lesion, and that in the background. In figure ITT] a few examples of 87 images are 
shown. Image (b) corresponds to the maximum 87 counting rate for the human scanner 
of about 1000 cps (cf. figure Ej) scanned for about 20 min. At contrast c = 8 the 1 cm 
lesion is clearly visible, however the one of 0.5 cm diameter is not. When contrast drops 
to c = 0.2 (c), the detectability of lesions deteriorates significantly. On the other hand 
the 1 cm lesion is still visible even if the count statistics is reduced by almost order of 
magnitude (d). 

Detectability of lesions clearly depends on their size and contrast. Let us examine 
this dependence to assess the limits of three-photon imaging. In order for a lesion to 
be distinguishable against the background its total number of points (counts) above the 
background level, reduced by its possible statistical fiuctuations, should be significantly 
higher than the average statistical fiuctuation of the number of points in a background 
region of the same size as the lesion. This can be expressed as 

VoUbC- 6 ^VoUbic+l) = e ^/v^b. (7) 
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Figure 11. Examples of images obtained from 87 events in simulation of the human 
scanner. The water filled phantom (a) consisted of 6 high activity spheres of diameters 
4, 3, 2, 1.5, 1 and 0.5 cm embedded in the 20 cm sphere with background activity, 
which was placed at the centre of the human scanner (cf. Table The images are 
mean filtered projections of 1 cm thick central slice, (b): background count density 
n(,=240 counts/cm^, contrast c=3, mean filter kernel (square averaging window) size: 
mm (c): 716=240 counts/cm^, c=0.2, filter 17 mm, (d): ni,=30 counts/cm"^, c=3, 
filter 9 mm. 



Vo is the original volume of the lesion without the partial volume effect. In the image, 
however, the blurring due to finite resolution of the scanner has to be taken into account; 
the blurred lesion volume is denoted by V. The parameters ^ and .^1 specify what 
"significantly higher" actually mean, in other words they control the confidence level 
of our assessment. Here we assume ^ = ^1 = 2, which corresponds roughly to 95% 
confidence level. For a spherical lesion equation ((Tj) takes the form 

y|^cao^-^iy^(c+l)ao3 = ^^(ao + (T)3 (8) 

where oq is the original diameter of the lesion and a is the FWHM spatial resolution of 
the scanner. Equation (jH)) is polynomial and can be solved numerically for oq. In the 
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ao 



\ 



17 



(9) 



The numerical solution of (jH)) for our model human scanner (cr = 1.5 cm) is plotted in 
figure 1121 




Figure 12. Minimum size of spherical lesion to be distinguishable against the 
background as a function of contrast and density of counts in the background for 
a scanner of spatial resolution 1.5 cm FWHM (human). 

The formula (jH)) works quite well for the simulated images. The calculated minimal 
detectable sizes of lesions for the images in figure ^2 are: (b) 6.2 mm, (c) 20.8 mm, 
(d) 10 mm. For the model small animal scanner (spatial resolution 0.32 cm FWHM) 
with 2 cm water phantom (section 14. 2j) of activity 50 MBq the 87 counting rate can 
reach over 50 kcps (see figure ITUj) . With the much smaller phantom the count density is 
about 1.5 • 10^ c/cw? in a 20 min. scan. Then sub-millimeter detectability is achieved 
at contrasts > 0.7. At c = 3 lesions as small as 0.6 mm should be visible. 

The formula (jHJ allows to fine tune the trade-off between the resolution and number 
of counts represented by the relation shown in figure |H1 Depending on the contrast of 
lesions, it may be of advantage to narrow the 87 window, losing some of the counts 
(decrease n?,), but improving the spatial resolution (decrease cr) to achieve the best 
detectability. 
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4.4- Attenuation 

In 37 imaging we unfortunately do not have the advantage of relatively easy attenuation 
correction as in conventional 27 PET, where the attenuation along any line of response 
is constant, independent of the site of annihilation and can be measured directly by 
a transmission scan. For each point in an object and each particular combination of 
energies and emission directions of the three photons the attenuation factor will be 
different. The problem of attenuation correction is therefore rather complex, similar to 
for example scatter correction in conventional PET or SPECT. Using Monte Carlo 
simulations and having the map of attenuation coefficients for the object one can 
compute a map of attenuation factors for each point of the object averaged over the free 
parameters of 87 emission. The results for our water phantoms are shown in figure ^1 



0.9 




I ' ' ' ' ' 1 

2 4 6 8 10 12 

r [cm] 

Figure 13. Average probability of emission from a sphere filled with water (diameters 
4 cm and 20 cm) without interaction for a 27 annihilation pair and the three photons 
from a 87 decay. 

In general attenuation for 87 photons is higher than for 27 because of lower energies 
and longer effective paths. It is also strongly nonuniform, being significantly higher for 
points deep inside the body. 

5. Outstanding questions and directions of future studies 

In this paper we concentrated on the image quality and performance of the new 87 PET 
modality, leaving for future studies the question of biological and chnical significance of 
the information extracted from the three-photon annihilations. Precise measurements 
of 87 yields in biological samples are necessary to determine its variability in a living 
organism, as well as sensitivity to parameters like the level of oxygen. Our results will 
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allow to assess whether the new information can be extracted by PET systems based 
for example on CdZnTe detectors. 

In our simulations we neglected the effect of finite positron range and non-zero 
residual momentum of the annihilating positron-electron pair. Both effects add extra 
blurring of magnitude similar to that for the two-photon case. That is in most cases 
negligible in comparison to the effects of detector size and energy resolution. 

In our estimation of counting rates we did not consider the contribution of activity 
outside the field of view, which is important in conventional clinical PET. It would 
add extra random 87 counts. However, our calculations show (section 14. 2j) that for 
the assumed scanner parameters the rate of random 87 events is at least 8 orders of 
magnitude lower than that of true 87. Extra activity contributing to randoms would 
mean a shift on the randoms curves in figures l9l and ITOl along the activity axis typically 
by a factor of 2-4. It means that the random 87 rate would still remain very small 
compared to the true 87 and would not pose a real problem. The above is true for pure 
positron emitters, which are most commonly used in PET. If there are nuclear gamma 
photons accompanying positron emissions, the rate of triple coincidences and false 87 
events may increase significantly, although still most of them could be rejected due to 
the energy conservation condition (jT]). 

The limitations of 87 imaging seem to be in the first place the low number of 
counts, especially with high attenuation (large patients, areas located deep in the 
body), and the spatial resolution limited mainly by the energy resolution of detectors. 
While one can not realistically do much to increase the counting rates achieved in 
our model scanners (section 14. 2|) . the spatial resolution leaves plenty of room for 
improvement. One way would be to rectify the energy resolution of the detectors (see 
section l4.1.ip . Although the technology of room temperature semiconductor detectors is 
making constant progress, we probably can not expect a significant further improvement 
in energy resolution compared to that assumed in our simulations. An alternative is to 
use e.g. HP-Ge detectors offering resolutions of the order of 0.8 %. However, the other 
properties, and in particular the need for cryogenic cooling and cost make the choice 
rather impractical at least for clinical PET. Their use for a dedicated small animal 
PET system is, however, not excluded (Philips et al 2002). A much more feasible 
way to improve the resolution of 87 imaging is through the reconstruction process. 
In this paper we used the simplest possible procedure by solving equations ((T)) for 
each registered 87 event and producing a "set of points" image. The easiest way to 
improve resolution is by making use of the energy dependence of the positioning error 
(figure Ej) as indicated in section 14.1.51 Another possibility is to account for the unique 
non-symmetric shape of the point spread function which can be calculated for each 
combination of detected 87 photons. It could be incorporated in the framework of a 
statistical iterative reconstruction method like ML-EM (Shepp and Vardi 1982), known 
to produce superior quality images, however, for the price of significant complexity and 
computation time. 
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6. Conclusions 

We have studied thoroughly the main characteristics of three-photon imaging in positron 
emission tomography. Clearly, it is one of the directions for future development of this 
rapidly expanding imaging technique. The main prerequisite for the new modality to 
be implemented is a scanner based on high energy resolution detectors, like CdZnTe. 
This is no longer a futuristic dream. Several groups and manufacturers' laboratories 
are working on such devices, and probably within the next few years we will see first 
prototypes in practical use. 

Our simulations show that for typical scanner configurations, with currently 
available CdZnTe detector properties, good quality 87 annihilation images can be 
obtained. They do not match those from conventional PET in terms of spatial resolution 
and statistics, however they may contain distinct new information for example about the 
oxygenation of tissues. It could be obtained alongside any routine scan, e.g. FDG PET, 
using photons that are currently wasted, so it is certainly worth further exploration. 
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